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When a mixture of two kinds of grains differing in size and shape is poured in a vertical two- 
dimensional cell, the mixture spontaneously stratifies in alternating layers of small and large grains, 
whenever the large grains are more faceted than the small grains. Otherwise, the mixture spon- 
taneously segregates in different regions of the cell when the large grains are more rounded than 
the small grains. We address the question of the origin of the instability mechanism leading to 
stratification using a recently proposed set of equations for surface flow of granular mixtures. We 
show that the stable solution of the system is a segregation solution due to size (large grains tend to 
segregate downhill near the substrate and small grains tend to segregate uphill) and shape (rounded 
grains tend to segregate downhill and more faceted grains tend to segregate uphill). As a result, the 
segregation solution of the system is realized for mixtures of large-rounded grains and small-cubic 
grains with the large-rounded grains segregating near the bottom of the pile. Stability analysis 
reveals the instability mechanism driving the system to stratification as a competition between size- 
segregation and shape-segregation taking place for mixtures of large-cubic grains and small-rounded 
grains. The large-cubic grains tend to size-segregate at the bottom of the pile, while at the same 
time, they tend to shape-segregate near the pouring point. Thus, the segregation solution becomes 
unstable, and the system evolves spontaneously to stratification. 



I. INTRODUCTION 



One of the unusual properties [gjj2| of granular mate- 
rials |3]-^| is the size-segregation of mixtures when they 
are exposed to external periodic perturbations such as vi- 
brations or rotations [j9fll5| . Size-segregation also occurs 
when a mixture of grains of different size is simply poured 
onto a heap |l6 - 21 : the large grains spontaneously seg- 
regate near the bottom of the heap, whereas the small 
grains segregate near the pouring point at the top of the 
heap. 

Recently, it was shown p2| that when a mixture com- 
posed of grains differing not only in size but also in shape 
is poured in a "granular Hele-Shaw cell" (two vertical 
slabs separated by a gap ~ 5 mm) , a spontaneous strati- 
fication is observed. Granular mixtures of small-rounded 
grains and large-cubic grains stratify in alternating lay- 
ers of small-rounded and large-cubic grains parallel to the 
surface of the pile when they are poured in the cell. 

According to the experiments , the control param- 
eter for stratification appears to be the difference of the 
repose angles of the pure species 

5 = e 22 -e u , (i) 

where 9n is the angle of repose of the small grains, and 
6*22 is the angle of repose of the large grains. The strati- 
fication experiments used a mixture of grains of dif- 
ferent shapes (small rounded or less faceted grains and 
large cubic or more faceted grains). The repose angle of 
the smaller pure species is then smaller than the repose 
angle of the large pure species — i.e., 5 > 0. On the other 
hand, strong segregation but not stratification occurred 



p2fl when 5 < 0, (corresponding to a mixture of small- 
cubic grains, and large-rounded grains) [p3| . We notice 
that the angle of repose of the pure species does not de- 
pend on the size of the grains, and it is a function of the 
shape of the grains: the rougher the grains the larger the 
angle of repose. 

To describe the case of a single-species sandpile 
in a two-dimensional geometry, Bouchaud, Cates, R. 
Prakash, and Edwards (BCRE) f|||27| developed a 
novel theoretical approach. They introduced two coarse- 
grained variables: the local height of the sandpile, and 
the local "thickness" of the layer of rolling grains, and a 
set of coupled equations to govern the flow of the rolling 
grains and their interaction with the sandpile. Recently, 
de Gennes |2^] applied the BCRE formalism to the case 
of granular flows in a thin rotating drum |29J3(J , and very 
recently Boutreux and de Gennes (BdG) |31| treated the 
case of granular flows made of two species of different an- 



gles of repose. Makse, Cizeau and Stanley (MCS) [|32|j33 



reproduced stratification and segregation using a discrete 
model and the continuum approach developed by BdG. 
They showed that a "kink" mechanism |32||33| describes 
the dynamics of stratification in agreement with experi- 
mental findings |^2| . 

In this paper, we address the question of the origin of 
the instability leading to stratification. We study an- 
alytically segregation and stratification as observed in 
p2[ when the two species have different size and dif- 
ferent shape (or, in general, different angle of repose) 
p4j . We use the continuum approach of BdG [|TJ , and 
MCS f{2|]3"3] ] , to calculate the steady-state solution of the 
equations of motion for surface flow of granular mixtures. 
This solution shows the complete size-segregation of the 
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mixture with the large grains being found at the bottom 
of the pile. We then study analytically the conditions un- 
der which the instability leading to stratification occurs. 
Stability analysis show that the steady-state solution is 
stable under perturbations only when S < 0, while the 
segregation solution is unstable (leading to stratification) 
when S > 0, in agreement with experiments p2|] . 

The stratification instability can be seen as follows. 
There are two segregation mechanisms acting when pour- 
ing a mixture of grains differing in size and shape in a 
cell: 



granular mixtures. In Section III we calcu late the steady- 
state solution of the problem. In Section IV we perform 
a stability analysis, and in Section |y| we discuss the in- 
stability mechanism for stratification and we propose a 
phase diagram for surface flows of granular mixtures in 
light of these results. 



II. THEORY FOR SURFACE FLOW OF 
GRANULAR MIXTURES 



• Size-segregation: large grains tend to segregate 
downhill near the bottom of the pile, and small 
grains tend to segregate uphill near the pouring 
point, since large grains roll down easier on top of 
small grains than small grains on top of large grains 
(Fig. ga). 

• Shape-segregation: cubic grains tend to segre- 
gate near top of the pile, and rounded grains seg- 
regate near the bottom, since rounded grains roll 
down easier than cubic grains (Fig. |l|b). 

Thus, when pouring a mixture of small-cubic grains 
and large-rounded grains the segregation of the mixture 
results, since the small-cubic grains size-segregate and 
shape-segregate near the top of the pile, and the large- 
rounded grains size and shape-segregate near the bottom. 
This situation gives rise to the steady-state solution of 
the system. 

On the other hand, the stratification process arises as 
a consequence of an instability mechanism. For mixtures 
of large-cubic grains and small-rounded grains there ex- 
ists a competition between size-segregation and shape- 
segregation. Large-cubic grains tend to size-segregate at 
the bottom of the pile, while at the same time, they tend 
to shape-segregate at the top of the pile. Thus, the seg- 
regation solution becomes unstable, and the instability 
drives the system spontaneously to stratification. 



The theoretical study of surface flows of granular ma- 
terials was triggered by the wor ks of BCRE HH and 
Mehta and collaborators 6{J. In a recent theoretical 
study for the case of a single-species sandpile BCRE 
f26| , P7f proposed two coupled variables to describe the 
dynamics of two-dimensional sandpile surfaces: the local 
angle of the sandpile 9(x,t) (or alternatively the height 
of the sandpile h(x,t)) which describes the static phase 
(i.e., the grains which belong to the pile), and the local 
thickness of the layer of rolling grains R(x, t) to describe 
the rolling phase (i.e., the grains that are not part of 
the pile but roll downwards on top of the static phase). 
BCRE also proposed a set of convective-diffusion equa- 
tion for the rolling grains, which was later simplified by 
de Gennes M: 



dR(x,t) 

at 



OR 

dx 



r(R, 



(2) 



where v is the downhill drift velocity of the grains along 
x, assumed to be constant in space and in time. The 
interaction term T takes into account the conversion of 
static grains into rolling grains, and vice versa. The sim- 
plest form of r is ||26 



(3) 



r(i?,6») = 7 (6{x,t) -Q r ) R(x,t). 



a) Size-segregation 



b) Shape-segregation 





rounded 



FIG. 1. Two segregation effects acting when the grains dif- 
fer in size and shape, a, Size-segregation: large grains tend 
to segregate at the bottom of the pile, b, Shape-segregation: 
rounded grains tend to segregate at the bottom of the pile. 

In the following we take up each of these results in 
turn. The paper is organized as follows: In Section |n| 
we present the theoretical formalism for surface flows of 



Here r denotes the angle of repose (the maximum an- 
gle below which a rolling grain is converted into a static 
grain ) . The rate 7 > has dimension of in- 

verse time, and v/j represents the length scale at which 
a rolling grain will interact significantly with a surface at 
an angle slightly above or below the angle of repose [gsj . 
For notational convenience we do not consider the dif- 
ference between the angle and the tangent of the angle, 
i.e. 



<**>-£■ 



(4) 



The equation for h(x,t) follows by conservation 
dh(x, t) 



Of 



= -T(R,( 



(5) 



2 








L x 



FIG. 2. Diagram showing the variables used to describe 
the granular flow of mixtures in a granular Hele-Shaw cell or 
two-dimensional silo. 

Recently, BdG have extended the BCRE formal- 
ism to the case of two species. This formalism considers 
the two local "equivalent thicknesses" of the species in 
the rolling phase R a (x,t) (i.e. the total thickness of the 
rolling phase multiplied by the local volume fraction of 
the a grains in the rolling phase at position x), with 
a = 1,2 respectively for small and large grains. The 
total thickness of the rolling phase is defined as 



R{x,t) = Ri(x,t) + R 2 {x,t). 



(6) 



The static phase is described by the height of the sandpilc 
h(x,t), and the volume fraction of static grains 4> a (x,t) 
of type a at the surface of the pile. Here x is the longi- 
tudinal coordinate, and the pouring point is assumed to 
be at x — 0, and we consider a silo or cell of lateral size 
L (see Fig. %. 

The equations of motion for the rolling species are BL] 



dR a (x, t) 



— Vn 



8R r 



la , 



dt dx 
and the equation for h(x, t) follows by conservation 



dh(x, t) 
dt 



-ri-r 2 



(7a) 



(7b) 



Here v a is the downhill convection velocity of species a 
along x. The interaction term T a takes into account the 
conversion of rolling grains of type a into static grains, 
and the amplification of static grains a by rolling grains 
of type a or (3, T a is defined through a 2 x 2 collision 
matrix M Q ^ 



T a (9, R ai (j) a ) = ^ M a pRp 

13=1 



(7c) 



The elements of the collision matrix M Q( a characterize the 
interaction of a rolling grain of type (5 with a surface of 
static grains of type a, and they are determined by the 
local angle 9(x,t), and the concentrations 4> a (x,t). The 
concentrations of static grains at the surface of the pile 
<fi a (x,t) are given by 



<f>a(x,t) 



dh 



and 



(7d) 



(7e) 



The canonical form of the collision matrix is defined 
by taking into account a set of binary collisions between 
a rolling and a static grains Hj3^] 



M : 



ai(0)0i 
xi(6)6 2 



h(6) x 2 (0)(f>i 
a 2 (6)<f> 2 - b 2 {9) 



(8) 



This definition involves a set of a priori unknown collision 
functions contributing to the rate processes: a a (6) is the 
contribution due to an amplification process, (i.e., when 
a static grain of type a is converted into a rolling grain 
due to a collision by a rolling grain of type a) , b a (9) is 
the contribution due to capture of a rolling grain of type 
a, (i.e., when a rolling grain of type a is converted into 
a static grain), and x a {9) is the contribution due to a 
cross-amplification process, (i.e., the amplification of a 
static grain of type (3 due to a collision by a rolling grain 
of type a). 

BdG used a "minimal" form of the collision matrix 
to calculate the steady-state solution in the geometry of 
a two-dimensional silo for the case of mixtures of grains 
differing only in angle of repose. This solution shows a 
complete segregation at the low edge of the silo, and near 
this point the concentrations of static grains show power- 
law behavior. In their model, they consider a constant 
cross-amplification term x a (9) =const, and they consider 
the angle of repose of each species to be independent on 
the surface composition of the pile. A generalization of 
the minimal model of pl[ to include not only different 
surface properties of the species but also different size of 
the species with small size ratios (d 2 /di < 1.4) is consid- 
ered in fl35|| . 

Here, we use the BdG equations to study segregation as 
well as stratification of mixtures of grains differing in size 
and shape in a two-dimensional silo (Fig. ||). As in MCS 
p2"t , we focus on the dependence of the repose angle of 
every species on the composition of the surface 4>q{x, t). 
We use the following definitions of collision functions [p2| : 



a a {9) 
b a {9) 



where 



U[x] 



laa Il[9(x,t) -9 a (<f>p)} 
l aa U[9 a (tt> fi ) -9{x,t)\ 
7/3q n[9(x,t) - 9 (<t> )}, 



if x < 

x if x > 



(9a) 



(9b) 
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Here, the rates are j aa > 0, and the generalized angle 
of repose 9 a (4>p) of a a type of rolling grain is a continu- 
ous function of the composition of the surface <pp ]32] (see 
Fig. ||, where we also define a p as a ((j>p) for <pp = 1): 



Ufa) 

Ufa) 



Tufa 
m0 2 



711 

hi 



-mq>i + W22, 



(10) 



where m = 
difference 



2) 



02i- We have assumed the 



(11) 



to be independent of the concentration 4>2- We notice 
that the fact that grains 1 are smaller than grains 2 im- 
plies 9\{(j)2) > #2(02) for any given concentration <p2 (he. 
the small grains are always the first to be captured). The 
angular difference ip is thus determined by the difference 
in size between the grains; the larger the size difference 
the larger if>. Moreover, the fact that the grains have dif- 
ferent shapes implies 9n ^ 822, (for instance #11 < #22 
when the grains 1 are more rounded than the grains 2). 




region A C <j> " §<-Q 

' e 22 \ 



FIG. 3. Dependence of the generalized angle of repose for 
the two types of rolling grains on the concentration of the 
surface of large grains <f>% for a, S > 0, and b, S < 0, where S 
is defined in Eq. (jl|). We define 9 a/3 = B a {4>0 = 1) H- 

The collision matrix includes the following processes 
@ (Fig. 1). 

(a) Capture, b a {9): rolling grains are captured if the 
local angle of the sandpile 9(x, f) is smaller than the gen- 
eralized repose angle 9 a ((j)p). The capture is proportional 
to R a go). 



(b) Amplification, a a (9): if the local angle 9{x,t) is 
larger than the generalized repose angle 9 a ((f>fi), then 
some static grains of type a will be converted into rolling 
grains due a collision by rolling grains of type a. The am- 
plification rate is proportional to the concentration <f) a in 
the sandpile, and to R a . 

(c) Cross-amplification, xp(9): when static grains of 
type a are amplified by rolling grains of type (5. This 
cross- amplification occurs when the local angle 9(x,f) is 
larger than 9r{c\>q). 

Equation ([7aJ) now reads 



0R o 
dt 



dR a 

' dx 



b 2 (9) 



U[9{x,t) - a (<f>p)] <t> a R a 

- l aa H[9 a ((f>p) - 8(x,t)} R a 
+ U[0(x,t) - 9 (3 ((j> )} 4> a R 

(12) 



a 2 (9) X 1 




m 2 ) 



m 2 ) 



FIG. 4. Schematic plot of capture b a (0), amplification 
a a (6), and cross-amplification x a (6) functions used in the def- 
inition of the present model. These functions are expected to 
be continuous in a region near the angle of repose (as shown 
by the thin curves). However, when the ratio between the size 
of the grains is not close to one, then tp = #1(^2) — #2(^2) is 
large enough so we can approximate these functions by the 
forms shown in this figure. 



III. STEADY-STATE SOLUTION. 

We now calculate the steady-state solution of the equa- 
tions of motion for the two-species sandpile including the 
terms corresponding to cross-amplification M Q ^, a ^ f3, 
which were not considered in |32] | . 

We consider the geometry of a silo of lateral size L (Fig. 
||). We assume that the difference tp = 9\{4>2) — ^2(^2) is 
independent of the concentration (f>2, then ip — 9\\—92\ = 
9\2 — 922 (see Fig. ||). We set v\ = V2 = v, and 
7n = 722 = 1, and we keep 712 7^ 721- We seek a solu- 
tion where the profiles of the sandpile and of the rolling 
grains do not change in time. Since stratification is an 
oscillatory solution, stratification cannot be observed for 
the steady-state solution. We set 



dh 
~dt 



R° 



(13a) 
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and 



dR a (x) 
dt 



0. 



(13b) 



with boundary conditions R a (x = 0) = i?", and R a (x = 
L) = 0, and R° = i?? + R%. 

We first calculate the profile of the total rolling species 
R(x) = Rx(x) + R 2 {x). From ©, ©, and (|l|) we 
obtain 







dR(x) 
dt 



dR(x) R° 
ox L 



(14) 



so that the profile of total rolling species decays linearly 
with a;, 



R° 

Lj 



(15) 



From ([7c]), Jl2]), and ( 13a ) we obtain the concentra- 
tions 



0i(x) 



7 n[0i(0 2 )-0(aO] R x (x) 



vR°/L + jn.[0 - 0i]Ri(x) + 72in[6» - 2 }R2{x) ' 

(16a) 

TjlN^HWj ga(g) 

vi?°/L + 7 n[6> - 2 ]i? 2 (aO + 7 i 2 n[6» - 1 }R 1 (x) ' 

(16b) 



The equations for the rolling species are obtained from 
(H), @, and d||), 



dR a (x) 
dx 



-(vR°/L) <j> a (x), 



(17) 



Next, we divide the calculations in two regions: Re- 
gion A, where 2 (<j> 2 ) < < 0i(0 2 ), and Region B, where 
< 02(<fe) < 0i (02) (see Fig. |b). The steady-state 
solution of Eqs. (0) and @ shows a strong segregation 
pattern. 

Region A. At the upper part of the pile we find that 
only small grains are present (for < x < x m , with 
x m = R®L/R° — v/^ipYsee below). 

If 2 < < 0i, from © we obtain §l| 



01 0*0 = 1, 



h(x) = o. 



(18) 



Using ( p7| ) and ( |l8| ) we find the profiles of the rolling 
species 



R° 

R 1 (x) = R%-—x, R 2 (x) = R%. 

Lj 



(19) 



The profi le of the sandpile is obtained from (7b), (12) 
and (13a) 



vR K) 
L 



so that 



= 7(00*0 -0ix)Rx{x)+^x{0{x) -02x)Ri{x), 



(20) 



0(x)-0 n = Q . (21) 

L(i?y + R-J,"f 21 /"f)/R° - x 

This solution is valid when 0(x) > 0i(0 2 = 1) = 2 i. 
Then it is valid for x < x m , where x m is 



_!% T _v_ 

Xm ~ BP 7^' 



(22) 



Region B. At the lower part of the pile (x m < x < L), 
we find that, after a small region of the order of v/("fijj), 
mainly large grains are present. 



If < 2 < 0i, from (|7bJ), (|2|), and ( |l3a|) we obtain 



^ = [0{x) - 1 (cl> 2 )}R 1 (x) + [0(x) - 2 (<f> 2 )}R 2 (x), 

(23) 



and therefore 



0(x)- 2 (h) = 



-vR Q /{ 1 L)+j ! R 1 (x) 
R(x) 



(24) 



Inserting ( |24| ) in (16) we obtain the concentrations as a 
function of the rolling species 



R(x) 



vR° 



R 2 (x) ( jipL 

Mx) = -57-r ( i - -^pRi(x) 



R(x) 



(25b) 



We obtain the equations for the rolling species using 
Eqs. @ and (H), 



dRi(x) _ Ri{x) fR° R 2 ( 



dx 



R(x) V L 



dR 2 {x) R 2 (x) (R° R x (x 



dx 



R{x) V L 



(26a) 



(26b) 



where r = v/(p/ip). Setting u — R\/R we obtain from 
~26a" 



ru 1 = — (1 — «)«, 



and the solution is 
Ri{x) 



R(x) 1 + C exp[(x - x m )/r] 



(27) 



(28) 



where C is an integration constant. By considering the 
continuity at x = x m , we obtain C — R?L/(R° r) M2| . 
The profile of the pile is obtained from ( p4[ ) and ( J10| ) 



5 



00) - 2 



-m<pi(x) 



ipR^L/R - v/-y 
(L-x) ■ 



(29) 



We notice that the parameter r = v/(pjip) is expected 
to be of the order of the size of the grains, so for large 
system size r <C L, and C S> 1. Then the steady- 
state solution can be simplified, and we arrive to the 
following simpler forms when also considering that, as in 
the experiment, an equal volume mixture is used (i.e., 
R\ = R Q 2 =R°/2). 

Region A. Valid for < x < x m = L/2 — vj (jtp). 



i? l( z) = ii°(i-£) 



(x)-> 



'11 = 



(t> 2 {x) = i) (30a) 

R 2 (x) = ^ (30b) 
-v/j - -072i£/(27) 



(L/2)(l +721/7) -z' 



(30c) 
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FIG. 5. Steady-state solution ©-© for the two-species 
granular flow in a silo geometry. We use the following typi- 
cal experimental values for the phenomenological constants 
appearing in the theory L — 30 cm, tan(9n = 0.6, 

tan#22 = 0.5, twaip = 0.2, tanm = 0.1, tan<5 = —0.1, 
R° — 0.25 cm, 11 = 10 cm/sec, 7 = 20/sec, and 721 = 10/sec. 
a, Profiles of the rolling species, b, Profile of the concentra- 
tions, c, Profile of the sandpile. 



Figure g shows the profiles of the steady state solu- 
tion for typical experimental values. We use a cell of 
size L = 30 cm. In |2f| the values of the different phe- 
nomenological coefficients appearing in the theory were 
measured for a typical equal-volume mixture consisting 
of quasi-spherical glass beads of mean diameter 0.27 mm, 
and cubic-shaped sugar grains of typical size 0.8 mm. 



The mean value of the velocity of the grains falling 
down the slope is of the order v ~ 10 cm/sec. The 
rate 7 was also estimated in |25j to be of the order of 
7 ~ 20/sec. A typical value of the thickness of the layer 
of rolling grains is R° ~ 0.25 cm jl3). The difference 
tan?/; = tan 0ii — tan$ 2 i is of the order of 0.1 — 0.3 
The concentrations (Fig. ||b) show the strong segregation 
of the mixture; the mixing of the species is concentrated 



G 



only in a small region, of the order of v/(j tan ip) ~ 1.5 — 5 
cm, in the center of the pile. The top part of the pile is 
made of small grains so that the angle is approximately 
equal to On (Fig. ||c). Towards the center of the pile the 
angle decreases and it is equal to 9 2 \ at x = x m (#21 is 
the angle at which large grains start to be captured on 
top of small grains). Then the angle increases gradually 
to 6*22 towards the lower part of the pile made of large 
grains. The profile of R(x) (Fig. |5|a) behaves linearly 
with x, which is a result of the conservation of number 
of grains Eq. (13a). Same dependence is observed in the 
e-species sandpile [E8| . The exponential 



case of a sing! 
behavior of the concentrations and rolling species R a (x) 
near the center of the pile is expected, since we are solv- 
ing equations of the type dR a /dx ~ — (■yip/v)R a (x). We 
notice that the segregation solution (pOh-(pl[) is deter- 
mined by the fact that 0i((f> 2 ) > 9 2 (4> 2 ), i-e. by the fact 
that the small grains are trapped first since they have a 
larger generalized angle of repose for a given concentra- 
tion </>2 . Then, this solution shows the size-segregation of 
the mixture. 



the center of the pile {x ~ x m ), without changing the 
angle of the pile. We then analyze the short-time evolu- 
tion of this perturbation to the steady-state profile. The 
dynamical evolution of the additional large rolling grains 
are now governed by the following equation: 



dR 2 {x 1 t) 
dt 



0R 2 
dx 



+ 71 



22 )R 2l 



(32) 



where the repose angle of the large rolling grains is 
since the surface is made only of large grains. 



'22 



(a) 5>0 



initial segregation 



(b) 5<0 
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IV. STABILITY ANALYSIS 



In the experiments of P2| , it was found that there is a 
transient regime of segregation before the layers appear. 
The initial segregation regime turns into stratification 
only for mixtures with 8 > (Eq. (|l|)). At the onset 
of the instability leading to stratification, it is observed 
that a small amount of large grains is captured on top 
of the region of small grains near the center of the pile 
(Region A, see Fig. ||a). This leads to the appearance 
of the first layer of larger grains and then to the oscilla- 
tions characteristic of stratification. On the other hand, 
if #11 > 6*22 (8 < 0) the segregation profile remains stable 
(Fig. ^b) . This picture was also confirmed by the models 
proposed in Q . 

The first sublayer of large grains corresponds to the ap- 
pearance of the first "kink" ||^,|3j| . The kink is an uphill 
wave at which the rolling grains are stopped. This first 
"kink" is formed only of large grains (as was observed 
in |^2| and |32|). This kink appears to move uphill and 
creates an incipient layer of large grains. Then the newly 
arriving grains roll down to the bottom of the pile where 
a new kink is developed. At this point, the kink (formed 
by large and small grains) starts to move upwards and 
all grains are stopped at the kink, but the small-rounded 
grains are stopped first so that the result is a pair of 
layers with the small-rounded grains being found under- 
neath the large-cubic grains. 

Next, we analyze the stability of the steady-state solu- 
tion under perturbations, by assuming that the steady- 
state solution obtained in Section [II is valid for the ini- 
tial transient regime of the evolution. We perturb the 



steady-state profiles by considering that a small amount 
of large grains has been captured in the region A, near 



first layer of large 
grains appears on 
Region A 




segregation 
remains stable 



mixing 




oscillations start: 
stratification 




FIG. 6. Schematic process leading to a, stratification 
(S > 0), and b, segregation (8 < 0). 

We are interested in the behavior of the profiles near the 
center of the pile, where R\ is very small, so we can fo- 
cus only on the behavior of the large grains. We look 
for the short-time behavior of R 2 , so we can assume that 
the angle of the pile remains unchanged from its initial 
value. Then we replace 9(x,t) in ( p2[ ) by 9(x) given by 
the steady-state solution (|2l|), and we arrive to the fol- 
lowing equation: 



dR 2 {x,t) 
dt 



dR 2 
dx 



7<5 - 



where v' = i[>j2iLR 2 /R°, and I 
The solution of (|33|) is 



(£-x) 
= L(R\ 



R 2 (x,t) = {l-xf e^ 5t , 
with uj = (v + v')/v. 



R 2 (x,t), (33) 
^§72i/7)/i?°. 

(34) 
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According to the exponential factor e 



-7<St 



(34) 



Ri{x, t) decreases as a function of time when 5 > 0. This 
implies that more large grains are captured, and the per- 
turbation of large grains evolves into the first sublayer 
of large grains and then to stratification according to 
the experimental findings |^| . On the other hand, when 
8 < 0, R2{x,t) increases as a function of time, so that 
the large grains of the initial perturbation are amplified, 
the perturbation disappears, and the segregation profile 
remains stable. 

(a) 8 < 0; Grains I : small and cubic (empty) 
Grains II: large and rounded (solid) 



stable 
segregation 





(b) 8 > 0: Grains I : small and rounded (empty) 
Grains II: large and cubic (solid) 



cubic 




unstable 
segregation = 
stratification 




large 



rounded 



size-segregation 



shape-segregation 



FIG. 7. The instability mechanism for stratification is a 
result of two competing segregation effects. Size-segregation 
(first panels in a and b), and Shape-segregation (second pan- 
els in a and b). a, Stable segregation when the mixture 
is composed of small-cubic grains and large-rounded grains, 
and b, unstable segregation leading to stratification when the 
mixture is composed of small-rounded grains and large-cubic 
grains. 



V. STRATIFICATION INSTABILITY 



III 



shows 



The steady-state solution calculated in Sec. 
the size-segregation of the mixture, while the stability 
analysis of Sec. IV shows that the grains also shape- 
segregate according to the value of 8. The stable steady- 
state solution is achieved when the mixture is composed 
of small-cubic grains (Grains I in Fig. |^a) and large- 
rounded grains (Grains II in Fig. ^a). In this case, the 
segregation of the mixture results because size and shape- 
segregation act simultaneously to segregate the large- 
rounded grains at the bottom and the small-cubic grains 
at the top. On the other hand, when pouring a mix- 



ture of small-rounded grains (Grains I in Fig. ^Jd) and 
large-cubic grains (grains II in Fig. |7^>) an instability de- 
velops since the size-segregation and shape-segregation 
mechanisms tend to segregate the same grain in opposite 
regions of the cell. 

§=022 — d 2 /d 1 = 1, shape-segregation 




stratification 




strong segregation 



5=0, size-segregation 



FIG. 8. Phase diagram for two-dimensional flows of gran- 
ular mixtures, is the size of the large grains, di is the size 
of the small grains, #22 is the repose angle of the large grains, 
and #11 is the repose angle of the small grains. We distinguish 
the regions of stratification and strong segregation (as shown 
by the steady-state solution presented here). Weak segrega- 
tion is expected in the shaded region when the size ratio is 
close to one. This case is discussed in [tS5||. 



Qualitatively, the onset of the instability can be seen 
as follows. If a small amount of large grains is captured 
near the center of the pile where the angle of the pile is 
8 ~ 6*11 (Point A in Fig. ^), then the repose angle for ad- 
ditional large rolling grains is 6*22- Thus, if 8 ~ B\\ < 822 
(Point B in Fig. |a), more large grains can be trapped 
(since the angle of the surface is smaller than the repose 
angle), leading to the first kink of large grains and then 
to stratification. On the other hand, the perturbation of 
large grains disappears when 8 < 0. Since 8 ~ 6\\ > 822 
(Point C in Fig. |b), no more large grains are captured, 
the fluctuation disappears, and the segregation profile re- 
mains stable. 

In light of these results we propose the phase diagram 
shown in Fig. |^. The model presented here predicts a 
region of stratification (8 > 0) , and a region of strong seg- 
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regation (S < 0) for size ratios not close to one. In Fig. g 
we show the results of numerical integration of the equa- 
tions of motion using the following parameters: stratifica- 
tion: 5 = 0.25, -0 = 0.3, 7n = 722 = 1, 712 = 721 = 0.1, 
vi = V2 = 1, Ri = i?° = 0-5, and segregation: the same 
parameters except for S = —0.1. For other parameters 
such as vi ^ i>2, different y a p, or 612 - #11 7^ #22 - #21 
we obtain similar results. In addition, when the cross- 
amplification rates j a f3 are of the order of the amplifica- 
tion rates 7 QQ we also find oscillations in the center of 
the pile which decay exponentially when S < 0. 



VI. DISCUSSION 

In summary, we study analytically segregation and 
stratification in granular mixtures focusing on the insta- 
bility mechanism for stratification. We find the steady- 
state solution of the equations of motion which shows 
strong size-segregation with the small grains located at 
the top of the pile, and with the large grains located 
near the bottom of the pile. In the center of the pile we 
find a region of mixing where the concentration profiles 
behave exponentially with a characteristic region of mix- 
ing given by v/f^ip). This region can be of the order of 
1.5 — 5 cm and is observed experimentally. We also find 
that the steady-state solution is stable under a pertur- 
bation involving large grains trapped at the center-top 
of the pile, when 5 < (corresponding to large grains 
less faceted than small grains). On the other hand, the 
steady-state solution is unstable under the same pertur- 
bation when S > (corresponding to large grains more 
faceted than small grains). The stratification instabil- 
ity is related to the fact that, when S > 0, there ap- 
pears a competition between size-segregation and shape- 
segregation: the large-cubic grains tend to size-segregate 
at the bottom of the pile, while at the same time, they 
tend to shape-segregate at the top of the pile. Thus, the 
segregation solution becomes unstable, and the system 
evolves spontaneously to stratification. 
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